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Abstract 

We study the ground state and the phase transitions of the bilayered spin- 
s' antiferromagnetic Heisenberg model using the Schwinger boson mean field 
theory. The interplane coupling initially stabilizes but eventually destroys the 
long-range antiferromagnetic order. The transition to the disordered state is 
continuous for small 5, and first order for large S. The latter is consistent 
with an argument based on the spin wave theory. The phase diagram and 
phase transitions in corresponding model in fractional dimensions are also 
discussed. 
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I. INTRODUCTION 



Recently there has been considerable interest in quantum spin liquids, which are magnetic 
systems without LRO at low temperature. While in general, the ground state of quantum 
spin systems lack true LRO in ID, the ground state of the 2D Heisenberg antiferromagnet 
exhibits Neel ordering even for S = 1/2, albeit with an sublattice magnetization that is 
considerably decreased from its classical value. Since spin is quantized, the spin value 
cannot be decreased beyond 1/2, hence the model does not have a spin liquid ground state. 
On the other hand, when two planes of antiferromagnetic spins are coupled together p-pT|, 
and if the interplane coupling is strong enough, the ground state is easily seen to be one of 
valence bond solid of interplane singlets (IVBS). Thus, there should be a transition from the 
LRO Neel state to a spin liquid state as the interplane coupling is increased. It has been 
suggested that the unusual magnetic properties of YBCO, with its basic unit of a pair of 
coupled CuO planes, may be due to its lying close to this quantum transition 0. 

It is of interest to study the nature of this quantum transition. Within a non-linear 



sigma-model (NLSM) description, Haldane |T2[ has pointed out that for a single plane of 
spins, topological Berry phase terms exist which differ between half-integer, odd integer, 
and even integer spins. One way to understand this is to consider the degeneracy of the 
valence bond solid states which maximize the number of resonating plaquettes in each case 
(4- fold, 2-fold, and non-degenerate respectively). On the other hand, the mapping of the 
two-plane system to the NLSM does not yield a topological term, which is consistent with 
the valence bond solid state for two planes with large interplane vs. intraplane coupling 
being zero- dimensional like and non-degenerate. Since the 2+1 D NLSM has only one phase 
transition which is second order, this suggests the same for the 2D quantum Heisenberg 
antiferromagnet at T = 0. However, the NLSM mapping assumes slow variation on the 
scale of lattice spacing, and so additional disordered phases and/or first order transition 
cannot be ruled out conclusively. 

In this paper we investigate the ground state of the 2D bilayered Heisenberg AF for gen- 
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eral S using the Schwinger boson mean field theory with no additional approximation . 
Our calculation complements previous calculations for S = 1/2 only and/or using additional 
approximations, as well as a calculation using the related Takahasi bosons approach Q-Q. 
Our results show that a first order transition is favored by decreasing quantum fluctuations, 
i.e.. increasing S. In particular, in agreement with those previous works, the transition for 
S = 1/2 is first order. The critical S separating first from second order transition equals 
0.35. While increasing interplane coupling J± eventually destroys LRO, ordering is stabi- 
lized by small J±, and the critical S for no LRO is shifted from the single plane value of 
0.2 to 0.13. A simple argument using spin wave theory helps to explain why first order 
transition occurs for large S. Since quantum fluctuations increases both with decreasing 5* 
and decreasing d, we also study the dimensionality dependence of the "bilayer" hypercubic 
system. For S = 1/2, we find that for d < 1.86 , the first order transition is replaced by 
second order one. En route, we also calculate the S vs. d phase diagram for J± = 0, i.e. the 
hypercubic Heisenberg AF. In addition to contradicting the NLSM description by having 
the possibility of a first order phase transition, whenever the transtion is continuous, the 
Schwinger boson MFT gives an additional phase transition between two disordered phases 
for all S, corresponding to a jump in the ratio of short-ranged intraplane to interplane cor- 
relations. However, since the Schwinger boson order parameter is not related directly to any 
physical symmetry breaking, it is likely this does not constitute a real phase transition but 
a sharp cross-over in behavior. 

II. BILAYERED ANTIFERROMAGNET IN 2D 

We begin with a quick review of Schwinger boson mean field theory [0] as applied to the 
translationally invariant nearest neighbor Heisenberg antiferromagnet on a bipartite lattice. 
The Hamiltonian is 

H = Y^ JijSi ■ Sj] Jij > 0, (ij) = n.n. 
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In the Schwinger boson representation, spin operators in each lattice site are replaced by 
spin 1/2 bosons as follows: 

The number of bosons at each lattice site is subject to the constraint: 

a 

which can be implemented by introducing a Lagrange multiplier on each site. The Hamil- 
toninan can now be written as 

^ = -2 E JijAljAj + ^NzJS^ + J2 Hblabia - 2-5), 

where Alj = ^Y.ablab]a bi^ = bii, bii = —bi-\ for sites on one sublattice and bia — bia 
for sites on the other sublattice. Physically, the product A\jAij acts as the valence bond 
(singlet) number operator for sites (i,j). In the mean field approximation, this product is 
decoupled by the Hartree-Fock decomposition. In addition, the exact local constraint is 
relaxed to one for the average: 

leading to the mean field Hamiltonian 

Hmf = Eo + XY, blk^ -2J2 JijAj{Alj + Aij), 

ia (ij) 

where we have taken Aij — (Aij) to be real. 

First consider the case that all the bonds are identical by symmetry, and assuming no 
spontaneous dimcrization, then all A^j must be the same Aij = A, . In this case Eq = 
jNzJS^ — 2XNS + JA'^Nz, where z is the coordiantion number. Hmf can be diagonalized 
by going to momentum space and perfroming the Bogoliubov transformation: 
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Hmf ^Eo-XN + ^u;k(4«k + PIPi. + 1), 

k 

where lo^ = [A^ - {JAz-fkf]^, 7k = i E^e*'*-'' = Ejli cosh/d. At T = 0, the energy should 
be minimized with respect to A and A, yielding the set of self-consistent equations: 

5 + --— V ^ 



2 27V^(^2_^2). 
1,^ .1 



k (/^^-7^)^ 

where we define n = X/ (JAz). An essential point of the theory is that a non-zero mean field 
amplitude A, which gives rise to boson hopping, indicates short-ranged antiferromagnetic 
order. Long-ranged order is achieved if the hopping amplitude is sufficiently large to give 
Bosc condensation. This occurs when these eqs. cannot be satisfied by having > 1, in 
which case /i — 1, and the k — term gives a finite contribution when converting the 
momentum sums into integrals: 

^1 1 r d'^k 1 

S+-^ms + 



2 2J-^{27rY^l-jl)h' 

1 d'^k -fl 

A^m, + - (1) 

2 J-n {2nf (1 - ^2^2 

It has been shown that the condensate density is also the sublattice magnetization. 
For the Heisenberg antiferromagnet on a square lattice, it was found that Bose condensation 
occurs for all S > Sc, where Sc = 0.2, with a gapless linear excitation spectrum characteristic 
of spin waves. For S < fji > I, and there is an energy gap for excitations. Thus, for all 
physical values of S, there is AFLRO. 

On the other hand, if two such planes are coupled together antiferromagnetically, and 
the interplane coupling is very large compared to intraplane coupling, the ground state is 
obviously a valence bond solid of interplane singlets, and the intraplane correlation length 
is zero. Thus, there must be at least one phase transition as the interplane coupling is 
increased. We now analyze this for general S using Schwinger boson MFT. 

The Hamiltonian in this case is 
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where sums over n.n. on the same plane and sums over n.n. on different planes. 

Since there is still translational invariance, the mean field Lagrange multiplier will be the 
same on all sites. However, the lack of symmetry between intraplane and interplane bonds 
means two mean field ampitudes must be introduced for the bond decoupling. Letting these 
be A and B respectively, and taking them to be both real, the mean field Hamiltonian is 
now 

Hmf = Eo + X - 2 E(4 + ^^j) 

ia (ij) 

where Eq = 2NJS^ - 2XNS + AJA'^N + + J^B^N. As before, we diagonalize 

Hmf through the Bogoliubov transformation, giving 



Hmf = Eq- XN + ^^ka(aL"ka + P^Pka + 1), 

ko- 

The excitation energies are given by 

d 

uJk,a = [A^ - {2 J A cos h + J^Baf] ^ . 

i=l 

where a = ±1. Minimizing Hmf with respect to A, A, and B gives the self-consistent 
equations: 



2 2JV ^ (1,2 - ^lJ^ 

A _ 1 Tk Q, , X)i=i cos ki , 

where Ty^a = {J2i=i cosfcj + Qa)/d and Q = J±B /(2JA). Note that fx, the excitation gap, 
must be greater than or equal 1 + Q/d. In particular, in the case of bose condensation, the 
value of /i is fixed to 1 + Q/d and hence the summations in Eq. (0) turn out to be a function 



of the parameter Q only. The magnetization is calculated by solving the self-consistent 
equations with the summations converted into intergrals: 



where /xq = 1 + Q/d. These equations hold so long as they give > 0, ie., Bose 
condensation, otherwise Eqs. @ should be used with /i also as an unknown parameter. In 
principle, we can solve for Q and or /x. In practice, the form of these equations allow us 
to avoid this by plugging in an arbitrary values of Q into the equations to find out or /x, 
and then A and B. The self-consistency is then reduced to using the values of A, B, and Q 
to determine the value of (3 = Jj_/ J. 

The behavior of Q as a function of (3 for S = 1/2 is shown in Fig. and is representative 
of all S. Discounting the trivial solutions Q = and Q = oo, corresponding to independent 
planes and IVBS respectively, there are Q ^ solutions indicating both intraplane and 
interplane correlations. For small (3, there is only one solution, with Q increasing from 
with j3. For (3 > 4(5* + 1/2)^, a second branch of solution, beginning at infinity appears. The 
two solutions merge at some larger value of (3, and beyond that, only the trivial solutions 
remain. The significance of these solutions can be understood if we consider the energy 
E{Q) obtained by minimizing the energy with respect to all other parameters except Q. 
Then, the non-trivial solutions are extrema of E{Q). Thus, for (3 < = 4(3 + 1/2)^, the 
solution of Q corresponds to a global minimum in E{Q), and describes the ground state. For 
f3 > Po, the upper branch corresponds to a local maximum while the lower branch remains 
a local minimum. The local maximum begins at Q = oo at /5o, and moves towards the 
local minimum with increasing (3. Eventually, the two extrema merge into a saddle point at 
f32- Beyond that, E{Q) is strictly decreasing with Q. By continuity, this means somewhere 
between /Sq and (32, E{oo) must cross from being greater than E{Qi) to less than it, where 




(3) 
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Qi is the lower branch solution. Thus, at this value the ground state jumps from the 2D 
correlated state described by Qi to the interplane VBS state. 

We can solve for the value of rris at the non-trivial solutions. Initially, increases with 
increasing Q, but eventually will decrease, vanishing at some Qc- For sufficiently large S, Qc 
will belong to the upper branch (maximum energy) solution. More importantly, in the case 
where Qc lies in the lower branch, its (3 value changes from less than to greater than (3i with 
increasing S. Thus, the transition from the LRO'ed state to disordered state is second order 
for small S, but becomes first order for larger S. In the former case, there is a subsequent 
transition from a disordered state with finite Q, hence with both interplane and intraplane 
short-ranged correlations, to the Q = oo state with only interplane correlations. Along with 
the jump in Q is a discontinous jump in the gap. It is tempting to associate this jump as 
a transition from some disordered state associate with a single plane to the non-degenerate 
IVBS. More likely this transition is probably an artifact of the Schwinger boson MFT, and 
indicates a relative sharp drop in the intraplane correlation length and a sharp rise in the 
gap. This is similar to the finite temperature MFT solution for a single plane, where A, hence 
short-ranged correlation, drops to zero above some finite temperatures In the latter case 
of first order transition in sublattice magnetization, the ground state jumps from one with 
LRO to the IVBS state. Since this latter state should be the correct ground state only in 
the f3 goes to infinity limit, we interpret this as the MFT way of showing a transition into 
a disordered state with a very short intraplane correlation length. The behavior of and 
the gap A as a function of (3 is shown in Fig. ^ for representative values of S. Fig. 2c shows 
an example of reentrance, where LRO first develops with increasing /?, but is subsequently 
destroyed when /3 gets too large. This occurs for S smaller than approximately 0.2, the 
MFT value of S below which the ground state has no LRO, but greater than approximately 
0.13, the minimum value of 5* for LRO at some (3. 

The phase diagram of S vs. [3 is shown in Fig. |. For S < 0.13, the ground state is 
always disordered. For 0.13 < S < 0.2, the system undergoes first a disorder-order and 
then a order- disorder continuous transition with increasing (3. For S > 0.2, there is LRO 
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for /3 = 0, and only the order-disorder transition remains. This transition is continuous 
until it terminates at a tricritical point at 5* ~ 0.35, /3 ~ 2.92, beyond which the continuous 
transition is preempted by a first order transition. Thus, for S > 0.35, there are values 
of (3 where the LRO'ed state is not the ground state, but is nevertheless metastable. The 
continuous transition phase boundary remains metastable until S ~ 0.4, beyond which the 
TJis = 0+ state moves into the upper branch and becomes unstable. In all cases of S where 
a disordered ground state with finite Q exists, a subsequent "first-order transition" occurs, 
with a discontinuous jump in Q and the gap A. As mentioned above, we interpret the jump 
as unphysical, and represents in reality a relatively sharp drop in the 2D correlation length. 

We can understand why large S favors a first order transition quite simply in terms of 
spin wave theory The Neel state energy is Ejy = S'^{2Jz + J_l) while the energy of 

the IVBS state is Ey = J_lS{S + 1). Equating the two implies an estimate for the first 
order transition at (3i of the order of S for large S. Within spin wave theory, the sublattice 
magnetization is given by rris in Eq. (Q) with B/A = 1. For large (3, the integral on the 
LHS scales as ^/J3. If we set rris = as an estimate for the critical value Pc for continuous 
transition, then /3c is of order S"^. Thus, for large S, (3i is much less than /3c. 

Within MET, the tricritical point and even the metastable continuous transition bound- 
ary occurs below the minimum physical value of S = 1/2. Thus, a first order transition is 
predicted for all physical systems described by the model. In fact, the sublattice magneti- 
zation jump at transition for 5* = 1/2 is about 30% of that at /3 = 0, clearly contradicting 
the results of numerical work on the model for 5* = 1/2, which supports a continuous tran- 
sition in the same universality class as the finite temperature transition of the 3D classical 
Heisenberg model. On the other hand, there is no reason to expect the Schwinger boson 
MET to give the exact answer, so the true position of the tricritical point might very well 
be above 5 = 1/2. This is particularly so since by relaxing the local constraint to a global 
one, unphysical states are included in the mean field solution, and the ME energy is not 
even variational. Thus, using these ME energies to find the position of first order transition 
is necessarily suspect. Nevertheless, we believe the prediction of larger S favoring a first 

9 



order transition to be correct, and the nature of phase transition in the bilayer system is 
non- universal. For example, the transition for S = 1/2 may become first order if there is 
a sufficiently large next nearest neighbor ferromagnetic interaction. Conversely, first order 
transition may become continuous if frustration is introduced. In other words, the value of 
S at the tricritical point can be changed by enlarging the parameter space. The seeming 
contradiction to the fact that the 2+1 D NLSM has only continuous transition is resolved 
by noting that the mapping of the Heisenberg model into the NLSM is legitimate only if the 
correlation length is long, which does not have to be the case of the disordered state close to 
a first order transition. Also of interest is that with sufficient frustration, the single- layered 
system can be disordered for S* = 1/2 or other physical values, and the reentrance behavior 
for small S discussed above in the bilayered system can be physically observed. 

Within our MFT and according to general arguments, first order transition implies the 
existence of metastable states with finite sublattice magnetization. This may lead to ob- 
serveable dynamics characteristic of macroscopic quantum tunneling. It would also be of 
significance with respect to Monte Carlo type numerical calculations due to problems 
of being "stuck" in the metastable minimum. For example, the first order transition may 
be missed if the metastability persists till the would-be continuous transition. 

III. EXTENSION TO THE FRACTIONAL DIMENSIONS 

We have seen that for the 2D bilayered square lattice antiferromagnet, Schwinger boson 
MFT shows the physically interesting case of 5 = 1/2 as undergoing a first order transition. 
Since the continuous transition is favored by small S, hence increasing quantum fiuctuations, 
one way of getting a continuous transition in MFT for S = 1/2 is to go to a dimension below 
2. In this section, we show that indeed this is the case. En route, we present the phase 
diagram of S vs. d (Fig. ^) for a single layer. These non-integer dimension results may 
be relevant to the physics of the Heisenberg antiferromagnet on percolating clusters, which 
have fractal dimensionality. 
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We first perform tlie Scliwinger boson MFT for a single hypercubic lattice in d dimension 
(square lattice for d = 2). The self-consistent equations (Eq. ([1|)) depend on d only through 
the momentum sum, which must be analytically continued to non-integer dimensions. This 
can be done by using the gaussian identity: 



In this form, the analytic continuation to arbitrary d is obvious (see Appendix). In Fig. 
the result for S* = 1/2 is shown. As d is decreased below 2, rris decreases and vanishes at 
some critical dimension dc = 1.46. Below dc, the excitation spectrum has a gap, which rises 
with decreasing d. These behaviors are representative for all S. However, it is known that 
for the simple Heisenberg Hamiltonian, 1/2 integer spin chains and integer spin chains are 
intrinsically different in that the former should be gapless, which can be understood as due 
to the presence of a topological term the appropriate NLSM. Thus MFT must break down 
for 1/2 integer spins even qualitatively as d gets sufficiently close to 1, and A must decrease 
again. 

Next we generalize our bilayer calculation to two coupled hypercubes in d dimension. 
The analytic continuation of Eqs (^ can again be done using the gaussian identity (^. As 
expected, a continuous transition can now be observed within MFT if d is reduced sufficiently 
from 2. For d < 1.46, the critical dimension for LRO for a single hypercube, the reentrance 
seen with increasing f3 for S < 0.2 is seen for = 1/2. So far, we have concentrated 
on lowering the dimension from 2. Of course, raising it would have the opposite effect. 
For example, for two 3D hypercubes, the 5 = 1/2 transition would be strongly first order 
within MFT, and so even taking into account inaccuracy of MFT, strongly implies a first 
order transition. 

Finally we discuss the critical phenomena of the continuous transition of this model. 
Analyzing Eqs. (2) and (3) close to the transition, we find the staggered magnetization 





to rewrite Eq. (|l[) as 




(5) 



11 



vanishes linearly, while the gap vanishes as {(3 — /3c)* with s = l/{d — 1) ioi d < 3, and 
s = l/2for(i>3 (there are logarithm corrections ai d = 3). These MF exponents are the 
same as those for the finite temperature transition of a single hypercube with the substitution 
d —>■ d + 1, refiecting the quantum nature of the present transition. 

IV. APPENDIX 

Notice that in Eq. (|^), the integral inside the bracket is a modified Bessel function of 
the first kind Io{x'^ — y"^). Therefore the first equation of Eq. (0) can be written as: 

S+l = m, + - r dxdy{e-^-'+y'^Io{x' - y'))". 

2 7C Jo 

Take the transformation: 

2 2 2,2 

u = X — y , V = X + y 

with V > \u\, for the integration variable will give us: 

1 2,d c~^'^ 

S+-=ms + — du{I^{u)Y / dv 

I TX J-oo j\u\ yv^ — u^ 

= rus — / du{Io{u)YKo{ud). 

IT JO 

It reduces the final formula into a single integral of Jq and Kq, modified Bessel function of 
the second kind. In this form, the integral indeed converges much faster than the original 
form in Eq. @ and consequently save much of the computation time. The same trick is 
also applied to both interplane and intraplane mean field equations. 
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FIGURES 

FIG. 1. The behavior of Q as a function of /? for S = 1/2. It is representative of all S. /3o = 4 
and P2 ~ 4.36. The star symbol stands for the location of 1st order transition in which f3 ~ 4.25. 

FIG. 2. Magnetization (solid line) nis and spin gap (dashed line) A//3 as a function of (3. (a) 
shows 1st order transition with the absent of 2nd order transition for 5 = 1/2 while both transitions 
are observed for smaller S in (b). (c), reentrance of magnetization occurs for S < 0.2. 

FIG. 3. Phase diagram of 5 vs /? for d = 2. 1st order transiton exists as long as (3 < 2.92 in 
which 1st order transition coincides with 2nd order transition (star symbol). The zero magnetiza- 
tion line (solid) terminates at a tricritical point, /3 « 3.38, beyond which the zero magnetization 
state is no longer stable. 

FIG. 4. Phase diagram of 5 vs d for a single layer. The curve goes to infinity as d tends to 1. 

FIG. 5. Magnetization iris and spin gap A for 1 < d < 2 and S = 1/2. The critical dimension 
that separates order and disorder state is 1.46. 
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